function cfdRunFalseTransientCase
%--------------------------------------------------------------------------
%
%  Written by the CFD Group @ AUB, Fall 2018
%  Contact us at: cfd@aub.edu.lb
%==========================================================================
% Routine Description:
%   This function cfdRuns steady state case
%--------------------------------------------------------------------------

cfdPrintHeader;

% Setup algebraic coefficients
theCoefficients = cfdSetupCoefficients;
cfdSetCoefficients(theCoefficients);

% Setup assembly fluxes
cfdSetupFluxes;

% Initialize cfdRun time
cfdInitTime;
cfdInitDirectories;

% Pre-updates (necessary on startup)
cfdUpdateFieldsForAllBoundaryPatches;% 处理边界上的值
% cfdUpdateGradients;% 计算梯度(但好像之后都没用到phi.grad的值？)
cfdUpdateScales;% 计算array里的最大最小值

% Start steady false transience loop
totalNumberOfIterations = 0;
while (cfdDoFalseTransientLoop)
    % Update number of iters
    totalNumberOfIterations = totalNumberOfIterations + 1;   
    
    % Time settings
    cfdUpdateRunTime;

    % Print
    cfdPrintIteration(totalNumberOfIterations);
    cfdPrintResidualsHeader;
    
    % Updates
    cfdUpdatePrevIter;% 将当前Region.fluid.var里的值赋给Region.fluid.var.prevIter
    % (作为上一迭代步的值存储下来),这里面还包含了变量的梯度等信息
    cfdUpdateProperties;% 根据当前计算结果，更新物性
    
    cfdUpdateGradients;% 但似乎挪过来之后结果没什么变化
    
    % Navier-Stokes and pressure correction (求解动量和连续方程)
    cfdAssembleAndCorrectNSSystem
    
    % Energy
    cfdAssembleAndCorrectEnergyEquation;
    
    % Post actions
    cfdPlotRes;
    cfdPostResults(totalNumberOfIterations);    
	cfdWriteResults(totalNumberOfIterations);
end